Ro- vibrational analysis of the XUV photodissociation of HeH+ ions 
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We investigate the dynamics of the photodissociation of the hydroheUum cation HeH+ by XUV 
radiation with the aim to establish a detailed comparison with a recent experimental work carried 
out at the FLASH free electron laser using both vibrationally hot and cold ions. As shown in previ- 
ous theoretical works, the comparison is hindered by the fact that the experimental ro-vibrational 
distribution of the ions is unknown. We determine this distribution using a dissociative charge 
transfer set-up and the same source conditions as in the FLASH experiment. Using a non-adiabatic 
time-dependent wave packet method, we calculate the partial photodissociation cross sections for 
the n = 1 — 3 coupled electronic states of HeH^. We find a good agreement with the experiment for 
the total cross section into the He -I- H dissociative channel. By performing an adiabatic calculation 
involving the n = 4 states, we then show that the experimental observation of the importance of 
the electronic states with n > 3 cannot be well explained theoretically, especially for cold (v = 0) 
ions. We also calculate the relative contributions to the cross section of the E and H states. The 
agreement with the experiment is excellent for the He^ -I- H channel, but only qualitative for the 
He -I- H^ channel. We discuss the factors that could explain the remaining discrepancies between 
theory and experiment. 

PACS numbers: 33.80.Gj 



1. INTRODUCTION 



The hydrohelium cation HeH+ has always attracted 
the attention of theorists. With only two electrons and 
two nuclei, it is one of the simplest heteronuclear molec- 
ular ions and it has been widely used as a benchmark 
for theoretical methods. In addition, being composed of 
the two most abundant species in the universe, it has 
many astrophysical applications. HeH+ is thought to be 
the first molecular ion to form in the early universe, by 
radiative association of He and H+ [l| . It has been pre- 
dicted to be abundant in various astrophysical objects 
such as planetary nebulae 0, helium-rich white dwarfs 
or metal-poor stars Q, although to this day it still 
eludes astrophysical detection. Several studies have been 
devoted to the search for IIeH+ in the planetary nebula 
NGC 7027 (see |^ for a review) , but it has appeared that 
its detection would be difficult. The abundance of IIeII+ 
in astrophysical environments is determined by the bal- 
ance between creation and destruction mechanisms. In 
the presence of a central star, such as in the case of plane- 
tary nebulae, the dominant destruction process of HeH^ 
is its photodissociation by energetic UV photons. It can 
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occur trough two dissociative pathways: 

HeH+(XiS+) + hv — > (HeH+)* — > He(lsnZ ^L) + H+ 

— > He+(ls) + H(7iO 
(1) 

As a consequence, several theoretical studies of the pho- 
todissociation process were realized. The first theoretical 
estimate of the photodissociation cross section is due to 
Saha et al. Q , who considered the transition to the first 
excited A ^E"*" state. The rate constant for the photodis- 
sociation process into the A state was first calculated by 
Roberge and Dalgarno using more accurate molecular 
data, while a calculation including the B state and 
the lowest ^11 state was done by Basu and Barua [s*]. 

In 2007, the first experimental data on the photodisso- 
ciation of IIeII+ in its ground X state were obtained 
using the free-electron laser (FEL) FLASH in Hamburg 
*9], operating in the XUV at a wavelength of 32 nm 
(corresponding to an energy of 38.74 eV). One of the 
main interests of free electron lasers is the possibility to 
operate at wavelengths below 100 nm, which allows to 
probe excited electronic states of molecules. The exper- 
iment made at FLASH consists in a crossed beam pho- 
todissociation of HeH+ with three-dimensional reaction 
imaging. The events leading to neutral helium fragments 
were detected, and the kinetic energy release (KER) of 
these fragments was analyzed, along with the angular 
orientation of the dissociating molecule with respect to 
the photon polarization. The photodissociation cross 
section was also determined. This experiment on the 
HeH+ system revealed the importance of n > 3 states in 
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the photodissociation, as well as the dominance of pho- 
todissociation perpendicular to the laser polarization and 
therefore the major role of the 11 states during dissocia- 
tion. This demonstrated the necessity of taking into ac- 
count the perpendicular orientation generally neglected 
or incompletely considered in previous theoretical mod- 
els. Following the publication of these results, two theo- 
retical studies investigating the role of excited states in 
the photodissociation process were realized. Dumitriu 
and Saenz calculated the cross section using a time- 
independent approach and neglecting the non-adiabatic 
couplings. Sodoga et al. [Ill computed the cross section 
using a time-dependent approach and showed the impor- 
tance of these couplings in the photodissociation process, 
but considered only states up to a principal quantum 
number n = 3. Despite these extensive calculations, an 
agreement with the experiment was difficult to obtain, 
one of the main reasons being that the experimental ro- 
vibrational distribution of the HeH+ ions was not known. 
While both theoretical works concluded that the cross 
section was not significantly affected by rotational exci- 
tation, the effect of vibrational excitation was found to 
be substantial. The most important unknown quantity 
is therefore the vibrational distribution of the ions in the 
experimental source. 

In this paper, we build on the work of Sodoga et al. 
with the aim to achieve a more detailed comparison with 
the experiment. By reproducing the ion source condi- 
tions, we determine the ro-vibrational distribution of the 
experimental ions in the FLASH experiment, and we use 
it to compare the theoretical and experimental cross sec- 
tions. In the meantime, Pedersen et al. published new 
results on the photodissociation of HeH+ [l^]. In this 
new study, the authors produced a source of cold ions 
by trapping the ions for about 100 ms before studying 
the photodissociation process. After such a long storage 
time, it can effectively be considered that all the ions are 
in the v = state. Moreover, the authors were able to 
detect events leading to the dissociation into He"'" + H. 
For hot (vibrationally excited) and cold {v — 0) ions, the 
authors measured the branching ratio between the two 
dissociation pathways in ([1} as well as the relative con- 
tribution of the and ^H states to the process. By 
investigating theoretically the photodissociation process 
for both hot and cold ions, we will show that there ex- 
ists significant discrepancies between the theoretical and 
experimental results. In particular, we will argue that 
states with n > 3 should not play a major role in the pho- 
todissociation process at the experimental energy. We 
will attempt to find the various factors that could ex- 
plain these conflicting results. 

This paper is organized as follows: in section [21 we 
quickly summarize the theoretical method used to treat 
the photodissociation process. In section [31 we describe 
the experimental method which allowed the determina- 
tion of the ro-vibrational distribution. In section [31 we 
compare our results with the experimental data, and we 
discuss the possible reasons for the discrepancy between 



theory and experiment. 



2. THEORETICAL FRAMEWORK 

We will only recall here the principal steps in the theo- 
retical description of the photodissociation process. The 
reader is referred to references [llj and |13| for a com- 
plete account of the theoretical methods used in this 
work. Our approach is based on the separation of the 
electronic and nuclear motions. The first step is to solve 
the electronic Schrodinger equation using ab initio meth- 
ods, giving access to the potential energy curves (PEC) 
and non-adiabatic coupling matrix elements. The pho- 
todissociation cross section is then obtained using a time- 
dependent approach, which consists in the propagation 
of wave packets on the coupled electronic states. 



2.1. Electronic structure of HeH"*" 

The ground state of the HeH+ molecular ion is a ^S"*" 
state. In the dipole approximation, the photodissociation 
can only occur through excited ^S"*" (parallel orientation 
of the field with respect with the molecular axis) and 
states (perpendicular orientation). 

We considered here all the n — 1 — 3 and "'^H 

states of the HeH+ molecular ion. The potential energy 
curves for these states have been calculated using the 
ab initio quantum chemistry package MOLPRO version 
2006.1 [iJl, as described in detail in [ll]- An adapted 
basis set consisting of the aug-cc-pv5Z basis set [l3| sup- 
plemented by one contracted Gaussian function per or- 
bital per atom up to n = 4 has been used. The PEC 
have been calculated at the state-averaged complete ac- 
tive space self-consistent field (CASSCF) level and are 
shown in figure [1] 

This approach allows us to compute the non-adiabatic 
radial couplings Fmm' , which are the matrix elements of 
the operator Or in the basis of the adiabatic electronic 
functions {Cm}: Fmm' = {Cm\dR\Cm,') ■ The non-adiabatic 
radial couplings are used to build the diabatic represen- 
tation [l^ , and the adiabatic-to-diabatic transformation 
matrix ©(i?) is the solution to the differential matrix 
equation 



dnPiR) + F( 



= 



(2) 



We have solved this equation starting from R — oo, where 
we have the initial condition D(cx)) = I. The diabatic 
potential energy curves are the diagonal elements of the 
matrix U'^ = D"! • U • D, where U is the matrix of 
in the adiabatic representation, while the off-diagonal el- 
ements of U'^ are the diabatic couplings. In the diabatic 
representation, the PEC cross and the couplings are a 
smooth function of the internuclear distance R. 

In our description of the photodissociation process, 
we also included the non-adiabatic rotational couplings. 
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FIG. 1: Illustration of the photodissociation of HeH"*". The 
ion is initially in a ro-vibrational level of the X ^'E'^ ground 
state characterized by an energy E^j and a wave function 
'4>vj{R) (in this example, t; = 0, J = 0). The photoabsorption 
is denoted by a vertical arrow, and the dotted line corresponds 
to the experimental energy {hv = 38.74 eV). We show the 
adiabatic potential energy curves of all the n = 1 — 3 and 
states through which the photodissociation can occur. 



which are the matrix elements of the operator Ly and 
connect states with A A — 1. Finally, we computed the 
dipole matrix elements as they govern the photodissoci- 
ation process (see below). 

In Ref. [T2II, the authors underline the importance of 
states with n > 3 during the photodissociation. We 
checked this assumption by introducing the calculated 
adiabatic PEC for the n = 4 and states ob- 

tained in Ref. [l^ in our calculations. However, as we 
were unable to compute the non-adiabatic couplings be- 
tween these states, we could only use these states to cal- 
culate the photodissociation cross section in the Born- 
Oppenheimer approximation (see section [ 



2.2. Computation of the cross section 

The initial wave packet is constructed on each of the 
excited states by multiplying the wave function of the 
initial state ipvjiR), with energy Eyj, by the matrix 
element /zio,mA' (R) of the dipole operator between the 
ground state (m = 1, A = 0) and an excited state m of 



lA'J' 
vj 



(3) 



The wave functions ipvj{R) and the energies E^j were ob- 
tained by solving numerically the vibrational Schrodinger 
equation in the potential of the ground state using a B- 
spline basis set method as in Ref. [T^ . 

The initial wave packet is propagated in time using the 
split operator [18|, which requires the knowledge of the 
diabatic representation. The propagator is composed of 
the potential and kinetic energy operators, the action of 
which is carried out in the coordinate and momentum 
representations, respectively. The switch between the 
two representations is done using the fast Fourier trans- 
form algorithm. The time step was 0.2 a.u., or 4.84-10~^ 
fs. We used a spatial numerical grid covering distances 
from 0.5 to 100 a.u., divided into 2^^ equally-spaced in- 
tervals. In order to avoid reflections of the wave packet 
at the end of the grid, we introduced an absorbing poten- 
tial in the asymptotic region at an internuclear distance 
Rc = 80 a.u. 

The total photodissociation cross section is obtained 
from the Fourier-transform of the autocorrelation func- 
tion C{t) = Em{^'vf'^'iR^O)\^Tf'^'iR,t)) fl^. How- 
ever, since the experiment distinguishes between the dis- 
sociating channels He -I- H+ or Hc+ + H, we will only 
look at the partial cross sections, i.e. the contribution 
from each of the excited electronic states. These are ob- 
tained by a method using the Fourier transform of the 
asymptotic wave packets on the electronic states at an 
internuclear distance i?cio, located in the asymptotic re- 
gion but before the absorbing potential. At each time, 
the wave packet is analyzed at the point R = R^o and the 
partial cross section as a function of the photon energy 
E = hv for each excited state is given by [20| 



{E) = 



-E\Al^,^'-^'\ 



(4) 



where /x is the reduced mass, km = y/2^{E — E^^) is 
the wave number in the electronic channel m with an 
asymptotic energy E^, u is the fine structure constant, 
ao is the Bohr radius, and 
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Starting from an initial level J in the ground state, 
the photoabsorption can excite the molecular ion into 
states with J' = J ± 1 or ^H states with J' = J ± 
0, 1. To obtain the photodissociation cross section for a 
given rotational level J, it is necessary to sum the various 
contributions as 



Sr^ SjA,J'A' 
J' 



2J' 
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.raK'.J'{E) (6) 



The Honl-London factors Sjjyjijyi indicate how the in- 
tensity of a transition is distributed among the rotational 
branches and are given in Ref. ^21i] . 



3. EXPERIMENTAL DETERMINATION OF 
THE RO- VIBRATIONAL POPULATION OF THE 
IONS 

In order to determine the vibrational and rotational 
distribution prevailing under the conditions of the exper- 
iment of Pedersen et al., we took advantage of the disso- 
ciative charge transfer set-up developed for the charac- 
terization of photoions [l^ . A keV beam of HeH+ ions is 
accelerated out of a duoplasmatron source operating with 
a mixture of hydrogen and helium at similar pressure and 
arc current as in [l2|. The beam passes through an ef- 
fusive potassium jet, where charge transfer takes place. 
As a result, excited HeH molecules are formed, that un- 
dergo predissociation towards the He -I- H ground state 
asymptote: 

HeH+ -I- K — > (HeH)* +K+ — ^ He -I- H -I- K+ (7) 

The fragments are detected in coincidence by a pair of 
position sensitive detectors (Quantar), giving access to 
the total kinetic energy release (KER). This process has 
been extensively studied by van der Zande et al. [23j , 
who have identified the contribution of three different 
molecular states to the predissociation yield, namely the 
A C and B ^H states. The mechanism is illus- 
trated schematically in figure [5] Due to energy conserva- 
tion, the kinetic energy release spectrum exhibits around 
8 eV a group of peaks corresponding to the ro- vibrational 
ladder of an unknown combination of the three states. 
In order to lift that major limitation of the method, we 
have set up an electrostatic ion beam trap consisting of 
a pair of mirror electrodes switched on and off at appro- 
priate times to allow for injection and extraction of an 
ion bunch. With typical beam lifetimes of about 40 ms, 
we were able to let the ions cool down to the vibrational 
ground state before they undergo the charge transfer pro- 
cess in a separate chamber located downstream. After 
vibrational cooling, the KER spectra are dominated by 
the predissociation of the v = level of the A, B and C 
molecular states, allowing for an accurate determination 
of their respective contribution to the signal, together 
with an effective rotational temperature. The extracted 
weights, i.e. 1 : 0.45 : 0.075 for the A, B and C state 
respectively, are then used to fit a rotational and vibra- 
tional HeH"*" population to the hot spectra, making use 
of the Franck-Condon matrices connecting the ionic and 
neutral states. The choice of a Boltzmann distribution 
of rotational states, convolved with the experimental en- 
ergy resolution, greatly simplifies the otherwise unstable 
fitting procedure. 

The vibrational population extracted from measure- 
ments performed with a duoplasmatron source opera ting 
under conditions similar to the FLASH experiment |12| 
is given in Table H] It is similar to the population ob- 
tained by van der Zande et al. [l^ who assumed an 
equal contribution of the A and B states, namely a dom- 
inant population oi v = Q with a monotonously decreas- 
ing population of higher states, limited to the first five 
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FIG. 2: Illustration of the decay of excited states of HeH 
following the charge transfer reaction (0 via photoemission 
or predissociation. The corresponding KER is shown on the 
right of the figure. 



vibrational levels. While this distribution is almost inde- 
pendent of the source parameters, an extensive rotational 
excitation is observed, which happens to be very sensitive 
to the type of ion source (ECR, duoplasmatron, hollow 
cathode) and its operating conditions, with temperatures 
ranging from 3000 to 5000 K. Such a hi gh t emperature, 



already reported by several authors [23|, l2j| , is not sur- 
prising since the HeH"*' creation process involves heavy 
particle collisions between vibrationally excited Hj ions 
and He atoms. Conversely, some vibrational quenching 
may occur in high pressure sources, although no clear 
connection could be established between source condi- 
tions and internal energy of the ions. The value quoted 
here for a duoplasmatron source, T = 3400 ± 300 K, falls 
in agreement with the temperature of 3100 K obtained 
by optical spectroscopy (23 |. 



4. RESULTS AND DISCUSSION 

4.1. Total cross section 

The total cross section for the photodissociation into 
He{lsnl ^L) + H+ starting from the initial state v — 
0,J — is shown in figure [21 together with the experi- 
mental value, a = (1.4 ± 0.7) • 10"^* cm^. We see that 
the calculated cross section is larger than the upper limit 
of the experimental error bars. This result is slightly dif- 
ferent than the one previously reported, the difference 
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TABLE I: Vibrational distribution of the ground X state 
of HeH"*" for the experimental conditions of Ref. and the- 
oretical cross sections for the two dissociating channels in ((T|) 
at hi/ = 38.74 eV. Error bars account for both the statisti- 
cal error and the systematics introduced by the simultaneous 
fitting of a rotational temperature and a vibrational distribu- 
tion. 
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being due to an error in the energy scale in Ref. 
and to the large variation of the cross section around the 
experimental energy (see figure [3]). 

In addition to the vibrational distribution of the ions, 
we present in table U the cross sections for dissociation 
into He ± H+ or He+ + H at the experimental energy. 
The cross section for the photodissociation into He ± 
H"*" weighted by the experimental vibrational distribu- 
tion is presented in figure S] At 38.74 eV, we obtain a 
cross section of cr = 1.78 • 10~^® cm^, close to the ex- 
perimental value. At the experimental energy, while the 
V — cross section is larger than the experimental value, 
it is balanced by the fact that the cross section is much 
smaller for all the other values of v. An interesting ob- 
servation is that at the experimental energy, the cross 
sections for i; > all have very similar values (ranging 
from 3.25 • IQ-^^ to 5.28 • IQ-^^ cm^), so that the main 
source of error in our theoretical estimation comes from 
the proportion of ions in w = 0. 

It should be noted that for v = or for the vibrational 
distribution of table HI we observe a dominance of the ^H 
states (perpendicular orientation) in the photodissocia- 
tion process. 



FIG. 3: Photodissociation cross section into He(lsn/ ^L) ± 
H+ for w = 0, J = 0. The contribution from ^2+ and 
states (parallel and perpendicular orientation, respectively) is 
shown, as well as the cross section averaged over an isotropic 
orientation of the field and the experimental data with error 
bars. 
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FIG. 4: Same as figure [3l but for the vibrational distribution 
of table Hand J = 0. 



4.2. Influence of the n > 4 states 

One of the experimental observations is the important 
contribution arising from states with n > 3 in the dis- 
sociation process, which was seen for both hot and cold 
ions. However, as indicated in Ref. [l^, the origin of the 
dominant role played by these Rydberg states remains 
unclear. 

In our calculations, we could only include states up 
to rt = 3, so that we cannot reproduce this experi- 
mental observation. This limitation of our approach is 
mainly due to the fact that we were unable to compute 
the non-adiabatic radial couplings between the n > 3 
states. However, as reported in [l^, we were able to 
obtain the adiabatic PEC of the n = 4 and ^H 

states, as well as their transition dipole moments with 



the ground state. With the PEC and dipole moments, 
we can calculate the photodissociation cross section in 
the Born-Oppenheimer approximation {i.e., neglecting 
all non-adiabatic couplings). While the partial cross sec- 
tion will be incorrect (since the n — 4 states are strongly 
coupled and that the couplings modify the final popula- 
tion of the electronic states), we can still study the effect 
of these states on the total cross section. 

The results are shown in figure[5]for ti = 0. We observe 
that the total cross section is modified by the inclusion of 
the n — 4: states, but that it is virtually unaltered at the 
experimental energy and below. This can be understood 
by looking at the PEC presented in figure [TJ when a 32 
nm photon is absorbed, the excitation occurs primarily 
in the n = 2 — 3 states. In figure [SI we show the rela- 
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FIG. 5: Total photodissociation cross section for v = 0, J = 
calculated in the Born-Oppenheimer approximation. Black 
line: using the n = 1 — 3 sates; red line: using the n = 1 — 4 
states. 
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FIG. 6: Relative contribution of the n = 2 — 4 manifolds to the 
total photodissociation cross section for n = 0, as a function 
of the photon energy. The n = 1 contribution is too small to 
be shown. 



tive contributions of the ^'E^ and states to the cross 
section, grouped according to the value of the principal 
quantum number n. Around the experimental energy, 
the principal contribution is due to the n = 3 ^11 states, 
followed by the n = 2 ^S"*" states. From these figures, it is 
difhcult to understand how states with n > 4 could play 
such a major part in the photodissociation process for 
cold ions. A possible explanation would involve a popu- 
lation transfer from the ti = 3 to the n > 3 states due to 
the non-adiabatic radial couplings. However, this expla- 
nation seems to be somewhat incompatible with the the- 
oretical results published on the charge transfer in ll{nl) 
+ He"*" (Is) collisions '25'], which showed little interac- 
tion between manifolds of different value of n. Moreover, 
by comparing the photodissociation cross section calcu- 
lated in the Born-Oppenheimer approximation or with 
the non-adiabatic couplings included, we observed only 
a small transfer of population between the n = 2 and 
n = 3 states. 

Another striking result of the experimental study is the 
fact that the importance of highly excited states seems 
to be independent of the initial vibrational excitation of 
the molecular ion and of the dissociating channel, as can 
be seen on figures 6 and 7 of Ref. [l^]- However, we ex- 
pect that when the value of v is increased, excited states 
with higher energy become accessible since the energy of 
the initial vibrational state, Eyj, is larger. In addition, 
the vibrational wave functions for high v cover a more 
important range of internuclear distances, so that the 
Franck-Condon region extends to higher value of R and 
more excited states become directly accessible (see figure 
[T]) . To illustrate this phenomenon, we show in figures [7] 
and [8] the contribution of the n = 4 states to the cross 
section starting from the initial state w = 1, calculated in 
the Born-Oppenheimer approximation. Compared to the 
V = case, we observe that at the experimental energy. 




Photon energy (eV) 

FIG. 7: Total photodissociation cross section for v — 1, J — 
calculated in the Born-Oppenheimer approximation. Black 
line: using the n = 1 — 3 sates; red line: using the n = 1 — 4 
states. 



the contribution of the n = 4 states is more important. 
In particular, the n = 4 ^H states account for about 25 
% of the total cross section. As a consequence, excited 
states with n > 4 should play a more important role in 
the photodissociation process for hot ions compared with 
cold ions, which is not observed experimentally. 



4.3. Branching ratios 

Pedersen et al. measured the branching ratio between 
the two dissociating channels in ([1]), tTjje++H/'''Ho+H+- 
They found a preference for dissociation into He+ + H 
for cold ions, and an almost equal contribution from both 
channels for hot ions. We calculated this branching ratio 
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FIG. 8: Relative contribution of the n — 2 — 4 manifolds to the 
total photodissociation cross section for i; = 1, as a function 
of the photon energy. The n = 1 contribution is too small to 
be shown. 
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Hot ions 


0.66 


0.96 ±0.11 


Cold ions 


0.42 


1.70 ±0.48 



TABLE II: Branching ratio (7He++H/''"He+H+ a-nd comparison 
with Ref. td. 



at the experimental energy for hot and cold ions by in- 
cluding the n = 1 — 3 states using the values presented in 
table m The results are presented in table |lll We observe 
an important disagreement between theory and experi- 
ment, particularly for cold ions, which cannot be easily 
rationalized. For hot ions, the cross section for dissocia- 
tion into the He ± channel is in good agreement with 
the experiment (see figure S]) , which would point to an 
underestimate of the He"*" ± H cross section. For cold 
ions, the experimental cross section is subject to a much 
larger uncertainty due to the smaller number of events 
detected. Unfortunately, the value of the cross section 
was not reported for these ions, so that a direct compar- 
ison with the experiment is not possible. Moreover, as 
can be seen in figures [3] and m the cross section for disso- 
ciation into He ± H"*" is very sensitive to a shift in energy 
around the experimental value. As there is an spread in 
the photon energy of ±0.6 eV in FLASH, we checked the 
dependence of the results presented in table against 
such a shift. For example, if the photon energy is low- 
ered by 0.6 eV, the ratio crjj(,+ _|_jj/o'Hc+H+ is increased to 
0.94 for cold ions and to 1.10 for hot ions. 

Pedersen et al. also estimated the contributions of 
and ^H states to the photodissociation process. Our re- 
sults are presented in table [HT] for hot ions. We observe 
an excellent agreement with the experimental values for 
states dissociating into He+ ± H, but a quantitative dis- 
agreement for the He + H+ channel. Qualitatively, we 



This work Experiment 



He ± H+ E 


10 % 


30 ± 2 % 


H 


90 % 


70 ± 2 % 


He+ ± H E 


39 % 


38 ± 3 % 


H 


61 % 


62 ± 1 % 



TABLE HI: Contribution of and ^H states to the pho- 
todissociation cross section for vibrationally hot ions. Com- 
parison with the experimental values of 1^ for J = 0. 



still confirm that the photodissociation in this channel is 
dominated by the H states. 

However, the results of table IIIII appear to contradict 
the above discussion. We found an good agreement with 
the experimental cross section for dissociation into the 
channel He ± H+, so that we would expect to obtain an 
agreement for the angular branching ratio as well. In- 
deed, the only mechanism allowing the transfer of pop- 
ulation from H to S states is via the non-adiabatic ro- 
tational couplings, but their effect was found to be neg- 
ligible. Moreover, to reproduce the experimental results 
of table In] for hot ions, we noted that it would be nec- 
essary to increase the cross section for dissociation into 
He+ ± H. However, the agreement between theory and 
experiment for the angular branching ratio is excellent. 
It would therefore be necessary to increase the cross sec- 
tion into He"*" ± H in an equal part for E and H states, 
which seems quite difficult to conceive. 

The angular branching ratios for cold ions are pre- 
sented in table lIVl As in the case of hot ions, the agree- 
ment is excellent for the He"*" ± H channel but it is only 
qualitative for the He ± H+ channel. In the latter case, 
we confirm the dominant contribution from the H states, 
but this contribution is found to be larger than in the 
experiment. 

Until now, we have assumed the effects of rotational 
excitation to be negligible, as found in Refs. [l3| and 
[ll| . However, due to the high rotational temperature of 
the ions (see section [3]), we decided to investigate its role 
in the photodissociation process for cold ions into detail. 
Even though the cross section does not depend strongly 
on J, an effect is expected for large J due to the change 
in the energy of the initial state. 

If we assume a Maxwell-Boltzmann distribution, the 
population pj{T) of the rotational level J with energy 
Ej is given by 

_ (2J±l)e-^.^A^ 
P^^^^ - E,,(2J'±l)e-^^'AT (8) 

where k is the Boltzmann constant. As indicated in sec- 
tion |31 the temperature is T ~ 3400 K. The cross section 
for a Maxwell-Boltzmann distribution of initial rotational 
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This work, J=0 This work, pj Exp. 



He + H+ S 


7% 


8 % 


24 ± 6 % 


n 


93 % 


92 % 


76 ± 6 % 


He+ + H E 


46 % 


50 % 


50 ± 3 % 


n 


54 % 


50 % 


50 ± 5 % 



TABLE IV: Contribution of and states to the pho- 
todissociation cross section for vibrationally cold (v — 0) ions 
and comparison with the experimental values of Ref. 
The results presented in the first column are for an initial ro- 
tational number J — 0, while those in the second column are 
for the Maxwell-Boltzmann distribution ((Sjl. 



states is then given as 

We computed the photodissociation cross sections aj for 
excited rotational states up to Jmax = 22. The angular 
branching ratios for the weighted cross section a are pre- 
sented in table IIVI We observe that the effect is indeed 
very small. Moreover, this is also the case for the total 
cross section and the ratio crjjj,+_,_jj/CTHc+H+- 

As J increases, the effect of the rotational couplings 
can also be expected to increase. This would modify the 
angular distribution, as the action of these couplings is to 
transfer population from S to 11 states (and vice versa). 
However, we observed that even for large values of J, the 
rotational couplings do not modify significantly the cross 
sections, so that the discrepancies between experiment 
and theory do not appear to be related to a mechanism 
involving these couplings. 

Finally, we checked the dependence of the angular 
branching ratio against a shift in energy of about 0.6 eV, 
corresponding to the uncertainty on the FLASH photon 
energy, as we noted (see above) that this could improve 
the ratio of the cross sections cjj(,+ ^fj/(THo+H+ • While the 
values presented in tables [IT] - IIVI are modified by such a 
shift, this does not resolve the discrepancies between our 
calculations and the experimental values. 

In Ref. [13, the experimental results are compared 
with the theoretical results of Dumitriu and Saenz [13] 
and an excellent agreement for the ratio (THe+-i-H/'''He-i-H+ 
is found for vibrationally cold ions. As pointed out by 
Pedersen et ai, in the work of Dumitriu and Saenz there 
is a clear separation in energy between states dissociat- 
ing into He + H+ or Hc+ + H (cf. figure 1 of Ref. [iQ:]). 
Therefore, these authors conclude that the non-adiabatic 
couplings will not affect the ratio between the fragments 
and only take into account a single non-adiabatic radial 
coupling between two n = 2 states. However, it was 
shown by Green et al. [1^ that the first part of the 
electronic spectrum up to n = 3 results from an alter- 
nation of states dissociating into He -I- H+ or He''" -I- 



H. This was confirmed and extended to n = 4 states 
in the recent study of Loreau et al. [T3| . Furthermore, 
the electronic states are strongly coupled, so that the 
non-adiabatic couplings modify substantially the partial 
photodissociation cross sections. 



5. SUMMARY AND CONCLUSIONS 

In this paper, we have expanded on a previous theo- 
retical investigation of the photodissociation of HeH+ in 
order to establish a detailed comparison with recent ex- 
periments for vibrational hot or cold ions. We cal- 
culated the photodissociation cross section by propagat- 
ing wave packets on the coupled potential energy curves 
of the excited states of the ion. In order to compare our 
theoretical results with the experimental data, we have 
measured the ro- vibrational distribution of the ions using 
the same source conditions as in the experiment. 

With the vibrational distribution, we have obtained a 
cross section for dissociation into the He -I- H+ channel 
close to the reported experimental value. To get a good 
agreement between theory and experiment, it was neces- 
sary to take the non-adiabatic radial and rotational cou- 
plings into account. On the other hand, our value for the 
cross section for dissociation into the other channel. He"*" 
+ H, is smaller than the experimental one. For cold ions, 
only the ratio of the cross sections for the two dissocia- 
tive channels, crfj(,+ ^jj/CTHc+H+ : has been reported. Our 
calculated value is in disagreement with the experimental 
estimate. To further compare the theory and experiment, 
it would be necessary to know the absolute value of the 
cross sections. For each dissociative channel, we have 
also calculated the angular branching ratios. For both 
hot and cold ions, the agreement is excellent for dissoci- 
ation into the He"*" -I- H channel but it is only qualitative 
for the He -I- H+ channel. This finding is difficult to ex- 
plain since the total cross section for He"*" -I- H is too 
small compared to the experimental value, while it is in 
good agreement for He -I- H+. The only way to modify 
the angular branching ratios while leaving the total cross 
section unchanged involves the non-adiabatic rotational 
couplings; however, their effect was found to be negligi- 
ble even for large J. Moreover, the effect of rotational 
excitation was found to be very small despite the high 
rotational temperature of the ions in the source. 

In our opinion, the most surprising conclusion of the 
experiment is the importance of states with a principal 
quantum number n > 3 both for hot and cold ions. While 
we were limited to the n = 1 — 3 states due to the dif- 
ficulty of calculating the non-adiabatic radial couplings 
for higher lying states, we could still compute the pho- 
todissociation cross section including the n = 4 states in 
the Born-Oppenheimer approximation (that is, neglect- 
ing the radial couplings). As shown in figures [5] and [6l the 
contribution of the n = 4 states is negligible for u = at 
the experimental energy. Moreover, we argued that the 
importance of the n > 3 states should increase with the 
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vibrational excitation of the ions, something that is not 
seen in the experiment. The excess kinetic energy associ- 
ated with higher lying initial vibrational levels (« 0.35 eV 
per level) and higher rotational temperature could eas- 
ily mask the downward shift of the dissociation energy 
with increasing principal quantum number of the prod- 
ucts. These counteracting effects require further model- 
ing, which lies beyond the scope of the present paper. 
This therefore remains an open question. 

Due to the complexity of the experiment, there is how- 
ever some uncertainty on the results which could affect 
the comparison. As noted by Pedersen et ai, direct ion- 
ization of HeH"*" into He"*" (Is) -|- H+ becomes possible 
at the experimental energy for vibrationally excited ions. 
This would lead to events appearing as H"*" , which would 
modify the branching ratios. However, this would not 
resolve the discrepancies observed for vibrationally cold 
ions. Another effect that would affect the determination 
of the branching ratios is the possible presence of the 
a state of HeH+ in the source. As it has a very long 
lifetime [13], this state (which lies about 11 eV above the 
ground state) would be present for both hot and cold 
ions. The absorption of a 38.7 eV photon by ions in 
this electronic state would lead to direct ionization, and 
thus to the observation of He+(ls) + H+. The a 
state has been previously observed in HeH+ sources [131 ■ 
Its presence was inferred from the low kinetic energy re- 
leased by dissociative charge transfer with an Ar target, 
and its disappearance with ion storage time. We have re- 
peated this experiment and replaced Ar with He, making 
the charge transfer process more resonant for the a '^S]+ 
state of HeH+ than for its ground state. No sizeable ef- 
fect was observed, which contradicts the attribution of 
the low KER to the a ^1]+ state. It seems thus unlikely 



that HeH+ is present in the a '^S+ in an amount that 
would significantly modify the above results. 

In conclusion, while we have been able to explain some 
of the experimental findings, there remains discrepancies 
between theory and experiment. Despite the simplicity 
of the system considered, the fact that a large number of 
excited electronic states play a part in the photodissoci- 
ation process makes an accurate comparison difficult to 
achieve. From an experimental perspective, it would be 
very interesting to measure the cross sections and branch- 
ing ratios at lower energy in order to make the com- 
parison with theory easier. In particular, we observed 
that the cross section changes very quickly in the region 
around the experimental photon energy. From a theo- 
retical point of view, the problem is complicated by the 
fact that the non-adiabatic radial couplings cannot be 
neglected in the treatment of the photodissociation. To 
investigate the role of the n > 3 states, it would therefore 
be necessary to calculate explicitly these couplings. 
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